arXiv:1503.02756vl [astro-ph.CO] 10 Mar 2015 


Draft version March 11, 2015 

Preprint typeset using DTfTjX style emulateapj v. 5/2/11 


ARE THE EFFECTS OF STRUCTURE FORMATION SEEN IN THE CENTRAL METALLICITY 

OF GALAXY CLUSTERS? 

Tamer Y. Elkholy 1,2 , Mark W. Bautz * 1 2 , Claude R. Canizares 1,2 

(Dated: March 11, 2015) 

Draft version March 11, 2015 

ABSTRACT 

A sample of 46 nearby clusters observed with Chandra is analyzed to produce radial density, tem¬ 
perature, entropy and metallicity profiles, as well as other morphological measurements. The entropy 
profiles are computed to larger radial extents than in previous Chandra cluster sample analyses. We 
find that the iron mass fraction measured in the inner 0.15i?5oo shows a larger dispersion across the 
sample of low-mass clusters, than it does for the sample of high-mass clusters. We interpret this 
finding as the result of the mixing of more haloes in large clusters than in small clusters, which leads 
to an averaging of the metal content in the large clusters, and thus less dispersion of metallicity for 
high-mass clusters. This interpretation lends support to the idea that the low-entropy, metal-rich gas 
of merging haloes reaches clusters’ centers, which explains observations of Core-Collapse Supernova 
products metallicity peaks, and which is seen in hydrodynamical simulations. The gas in these merg¬ 
ing haloes would have to reach the centers of clusters without mixing in the outer regions, in order 
to support our interpretation. On the other hand, metallicity dispersion does not change with mass 
in the outer regions of clusters, suggesting that most of the outer metals come from a source with a 
more uniform metallicity level, such as during pre-enrichment. We also measure a correlation between 
the metal content in low-mass clusters and the degree to which their Intra-Cluster Medium (ICM) is 
morphologically disturbed, as measured by centroid shift. This suggests an alternative interpretation 
of the large width of the metallicity distribution in low-mass clusters, whereby a metallicity boost in 
the center of low-mass clusters is induced as a transitional state, during mergers. 

Subject headings: galaxies: clusters: general, galaxies: clusters: intracluster medium, X-rays: galaxies: 
clusters 


1. INTRODUCTION 

In today’s Universe, galaxy clusters are the largest 
bound structures, ranging in mass from about 10 13 to 
10 15 Mq. Their formation is understood within the 
paradigm of hierarchical structure formation, where large 
structures form through the merger of smaller structures. 
Because of their large masses, the Intra-Cluster Medium 
(ICM) in galaxy clusters is heated and compressed, such 
that bremsstrahlung emission becomes an efficient radia¬ 
tion process. This X-ray emission from the ICM contains 
within it ample information about the chemistry and the 
dynamical state of the ICM, and by extension the clus¬ 
ter as a whole. In particular we focus in this study on 
two quantities, metallicity and entropy, because they en¬ 
code the integrated effect of various physical processes 
occurring in clusters. 

Most metals have been produced in stars, predomi¬ 
nantly through thermonuclear burning, or in processes 
resulting from the extreme conditions in supernovae 
(SNe). Therefore, studying metallicity in galaxy clus¬ 
ters probes processes that produce metals, such as star 
formation and SN rates, as well as processes that mix 
and distribute metals such as mergers and central AGN 
a ctivi ty. 

iRenzini et all ( 199 3) were the first to point out that 
one cannot account for all the observed metals in the 
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ICM assuming the current SN Type la rate per unit lu¬ 
minosity was the same in the past. This problem was 
studied with better SN data and more co mplex che mical 
evolut ion a nd ICM enrichment m odels bv lPortinari et all 
(120041) and iLoewensteinl (|200fil ). who both showed the 
need for more metal production in clusters, compared to 
what is observed in the fi e ld. B oth studies above in addi¬ 
tion to lNagashima et all (12005D . suggest that a different 
stellar Initial Mass Function (IMF) must operate in field 
galaxies than in clusters. For example, a top-heavy IMF 
might be necessary in galaxy clusters. 

Hydrodynamical models predic ting the meta llicity of 
the ICM have been reviewed by iBorgani et all (120081) . 
The results of these studies in terms of the effect of dif¬ 
ferent IMF’s on the level of enrichment vary, as some are 
able to match or even exceed the observed metallicity in 
clusters, while others cannot match it, depending on the 
many details and parameters of the models. 

Two classes of SNe are responsible for producing the 
metals in clusters. Core Collapse SNe (SNCC) result 
from massive stars, whose lives are short compared to 
the time since the peak of star formation, at a redshift 
z ~ 2 — 3. Type la SNe (SNIa) result from low-mass stars 
and thus can be delayed by billions of years from the time 
of star formation. If the above picture of a prompt metal 
injection from SNCC, followed by a more gradual injec¬ 
tion from SNIa holds, then we should expect this to af¬ 
fect the metallicity radial profiles. E ach class of SNe pro - 
duces a different set of elements (e.g. (Werner et al.ll2008l) . 
Both SN classes produce iron, while, for example, SNCC 
are the main contributors of oxygen, neon and magne- 
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sium. Since the rate of SNCC at any epoch is roughly 
proportional to the then-current rate of star formation, 
and since clusters form relatively late, we expect SNCC 
products to have enriched the proto-ICM before cluster 
collapse. The SNCC products are therefore expected to 
be homogeneously spread through the ICM. Conversely, 
SNIa products are expected to be nearer the center of the 
cluster, as the Brightest Cluster Galaxy (BCG) accumu¬ 
lates more and more stars wit h time. However, us i ng ver y 
deep Chandra observati o ns. iSand ers fc Fa bianl d2006h , 
Simionescu et al.1 (l2009afl . iSimionescu et al.l l)2010h and 
Million et al.1 (|20 1 il l measure a large amount of SNCC 
products in the core of nearby galaxy clusters, present¬ 
ing a challenge to the above picture. 

Because entropy remains unchanged under adiabatic 
processes, it measures any heat input or output to the 
gas. In the case of the ICM, the main process that 
changes its entropy is gravitational shock heating. Other 
processes affecting the ICM entropy include radiative 
cooling, heating by AGN and SNe, as well as processes 
such as turbu lent dissipation and conducti on. The reader 
is referred to lMcNamara fc Nulseni (2007) for a detailed 
review of the above processes. 

Gravity is expected to be the dominant force in clus¬ 
ter formation, and thus most of the entropy generation 
is expected to result from shock heating of infalling gas, 
as a cluster forms. This process converts gravitational 
potential energy into heat. Because S oc T/n 2 / 3 , we 
expect the characteristic entropy scale of a cluster to 
scale proportionally to its virial temperature, or as M 2 / 3 , 
where M is the cluster mass. Hydrodynami cal simula¬ 
tions (e.g. iVoit, et all 120051 iNagai et, alJl2007T) produce a 
normalization of the entropy expected from gravity, as 
well as a radial profile, which is proportion al to a power- 
law of radius, r 1 ' 1-12 . Observationally, iPonman et al.l 
(Il999h were among the first to measure the departure 
of the entropy from the expected self-similar scaling 
wi th cluster mass. The same result was l ater found 
in iLlovd-Davies et al.l (|2000l h IPonman et all (|2003T) . and 
iPratt et al.l (l20ldh 

Some studies have focused on cases where the ef¬ 
fects of both heating and metal enrichment can be de¬ 
tected. High-resolution observations of the centers of 
clusters find regions of enhanced metallicity expected to 
have been ejected from the central regions of cluster s. 
by the central engine fe.g. ISimionescu et all l2009al lbl; 
IKirkpatrick et all 120091: IQ Sullivan et all l20TTT) . For ex- 
ample. iKirkpatrick et al.l (l2011f) showed that the di¬ 
rection of elongation of an enhanced-metallicity region 
is correlated with the direction of cavities and radio 
emission axes, originating from the center . On larger 
scales, however, si mulations bv lBorgani et al.l (|2005l f and 
iShort et~ahl (1201311 suggest that winds from SNe have lit¬ 
tle effect on the overall entropy of clusters. There is, 
however, a shortage of attempts to test these results ob¬ 
servationally. 

The aim of this work is therefore to apply a systematic, 
spatially resolved study of metallicity and entropy to a 
large sample of galaxy clusters, to look for any possible 
relation between the metallicity and the entropy of the 
ICM. We also aim to produce a dataset of entropy pro¬ 
files (in the form of temperature and metallicity profiles) 
for the community to employ in various galaxy cluster 


studies, and make it available electronically. 

Throughout this work a ACDM cosmological model 
is assumed, where the Hubble constant is Hq = 72 
km/s/Mpc, the matter density parameter Ll m = 0.26, 
the dark energy density is Ha = 0.74, and the universal 
baryon fraction, fb = Llb/Ll m = 0.169, where Lib is the 
baryon density. 

2. DATA SAMPLE AND ANALYSIS 
2.1. Data Sample 

Our sample consists of bright clusters — present in 
both the HIFLUGCS and ACCEPT samples — which 
were observed with Chandra’s ACIS instrument out to 
at least O. 2 .R 500 , where R 500 is the radius enclosing an 
average density that is 500 times the critical density of 
the Universe at the redshift of the observed cluster. More 
precisely, we start with the extended HIFLUGCS sam¬ 
ple of galaxy clusters, which is a flux-limited sample of 
clusters with X-ray flux /A [0-1 — 2.4keU] > 2 x 10 -11 
erg/s/cm 2 (iReiprich fc Bohringeil 12002 1. We search for 
the HIFLUGCS clusters which have Chandra coveragcQ 
out to at least 0.2i? 5O o. Once the above sample is identi¬ 
fied, we search the ACCEPT data for the aforementioned 
HIFLUGCS subsample, where entropy is measured in at 
least 3 radial bins beyo nd O.2i? 50 o. The ACCEPT study 
dCavagnolo et al.ll2009lj measured the entropy profiles of 
all galaxy clusters observed by Chandra , up to August 
2008, and has made the data available onlhnfl Forty 
eight galaxy clusters satisfy the above selection criteria. 

For all spectral analysis, we discard all observations 
taken before January 29, 2000 as the Chandra focal plane 
temperature for that period was — 110°C or higher, which 
increases the level of background. Chandra’s calibration 
is in general better for later dates. Only one cluster, 
Abell 401, is excluded from our analysis because it was 
only observed prior to January 29, 2000. We also ex¬ 
clude Abell 2255 because its available observation was 
too short to yield enough photons. The complete analy¬ 
sis will thus be presented for a sample of 46 clusters. 

We show in Table [0 the observational details of our 
sample, and the Chandra observations used in this work. 
For each cluster, we show the label which will be used 
to designate it hereafter, in figures and tables. We show 
the cluster’s celestial coordinates, and its redshift. We 
list the Chandra Observation ID’s (OBSID) we use in this 
work. Some OBSID’s are shown in parentheses. These 
are the observations that were used in imaging analysis, 
and excluded from spectral analysis, as described above. 
Finally, we show the total exposure time used in the spec¬ 
tral analysis of each cluster. 

2.2. Data Preparation 

We use the Chandra Interactive Analysis of Observa¬ 
tions software, more commonly known as CIAC0, for 
analysis in this work. More precisely, we use Cl AO’s 
4.2 version. Data are reprocessed following the guide¬ 
lines in the CIAO analysis thread “Reprocessing Data 
to Create New Level=2 Event File”, using the tool 

3 For the purpose of selecting clusters, we use the values of -R 500 
measured in lReiprich fe Bohringeil (120020 . 

4 http://www.pa.msu.edu/astro/MC2/accept/ 

5 http: //cxc. harvard. edu/ciao/ 
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TABLE 1: Chandra Observations. 


Cluster 3, 

Label b 

RA C 

Dec d 


Obsids f 

T (ks) s 

Abell 119 

all9 

00:56:15.392 

-01:15:17.78 

0.044 

4180, 7918 

57 

Abell 1413 

al413 

11:55:17.986 

+23:24:17.82 

0.1427 

1661, 5002, 5003 

121 

Abell 1644 

al644 

12:57:11.772 

-17:24:33.68 

0.0474 

2206, 7922 

70 

Abell 1651 

al651 

12:59:22.188 

-04:11:45.80 

0.086 

4185 

10 

Abell 1689 

al689 

13:11:29.495 

-01:20:29.02 

0.184 

5004, 540. 6930, 7289 

181 

Abell 1736 

al736 

13:26:54.235 

-27:09:48.65 

0.0461 

10428, 10429, 10430, 10431, 

35 






10898, 10899, 10900, 10901, 







12026, 12028, 12029, 13106, 







13107, 13108, 13109, 13110, 







13111, 13112, 13113, 13412. 


Abell 1795 

al795 

13:48:52.668 

+26:35:30.73 

0.0616 

13413, 13414, 13415, 13416’ 

437 






13417, 5286, 5287, 5288, 







6159, 6160, 6161, 6162, 







6163, (494) 


Abell 1914 

al914 

14:26:01.072 

+37:49:32.97 

0.1712 

3593, (542) 

27 

Abell 2029 

a2029 

15:10:56.091 

+05:44:40.94 

0.0767 

10437, 4977, 6101, 891 

112 

Abell 2063 

a2063 

15:23:05.323 

+08:36:28.49 

0.0354 

5795, 6263 

27 

Abell 2065 

a2065 

15:22:29.060 

+27:42:34.33 

0.0721 

3182 

22 

Abell 2142 

a2142 

15:58:20.103 

+27:13:58.76 

0.0899 

5005, (1196, 1228) 

68 

Abell 2147 

a2147 

16:02:15.608 

+ 15:57:53.77 

0.0351 

3211 

18 

Abell 2163 

a2163 

16:15:46.519 

-06:08:50.57 

0.201 

1653, 2455, 545 

89 

Abell 2204 

a2204 

16:32:46.922 

+05:34:31.40 

0.1523 

499, 6104, 7940 

97 

Abell 2244 

a2244 

17:02:42.517 

+34:03:37.46 

0.097 

4179 

57 

Abell 2256 

a2256 

17:03:59.388 

+78:38:44.57 

0.0601 

2419, (1386, 965) 

35 

Abell 2319 

a2319 

19:21:09.997 

+43:57:18.82 

0.0564 

3231 

14 

Abell 2657 

a2657 

23:44:56.531 

+09:11:28.75 

0.0404 

4941 

16 

Abell 2734 

a2734 

00:11:21.616 

-28:51:17.98 

0.062 

5797 

20 

Abell 3112 

a3U2 

03:17:57.627 

-44:14:20.34 

0.075 

2216, 2516, 6972, 7323, 7324 

108 

Abell 3158 

a3158 

03:42:52.591 

-53:37:50.03 

0.059 

3201, 3712 

56 

Abell 3376 

a3376 

06:01:57.312 

-39:58:25.80 

0.0455 

3202, 3450 

64 

Abell 3391 

a3391 

06:26:20.780 

-53:41:32.98 

0.0531 

4943 

18 

Abell 3571 

a3571 

13:47:28.580 

-32:51:14.35 

0.0397 

4203 

34 






5751, 5752, 6292, 6295, 


Abell 3667 

a3667 

20:12:36.316 

-56:50:40.74 

0.056 

6296, 889, '(513) 

430 

Abell 3822 

a3822 

21:54:06.292 

-57:51:41.06 

0.076 

8269 

8 

Abell 3827 

a3827 

22:01:53.279 

-59:56:45.99 

0.098 

7920 

46 

Abell 3921 

a3921 

22:49:57.845 

-64:25:44.13 

0.0936 

4973 

29 

Abell 399 

a399 

02:57:51.557 

+ 13:02:32.43 

0.0715 

3230 

49 

Abell 400 

a400 

02:57:41.119 

+06:01:20.03 

0.024 

4181 

21 

Abell 4038 

a4038 

23:47:43.200 

-28:08:38.30 

0.0283 

4188, 4992 

40 

Abell 4059 

a4059 

23:57:00.933 

-34:45:34.44 

0.046 

5785, 897 

110 

Abell 478 

a478 

04:13:25.199 

+ 10:27:53.90 

0.09 

1669, 6102 

52 

Abell 539 

a539 

05:16:36.680 

+06:26:34.63 

0.0288 

5808, 7209 

43 

Abell 644 

a644 

08:17:25.392 

-07:30:48.38 

0.0704 

10420, 10421, 10422, 10423, 
2211 „„„„ 

49 

Abell 754 

a754 

09:09:21.084 

-09:41:05.78 

0.0528 

10743, 6793, 6794, 6796, 

6797, 6799, (577) 

187 

Abell S 405 

as405 

03:51:29.787 

-82:13:21.26 

0.0613 

8272 

8 

Hydra A 

hyda 

09:18:05.876 

-12:05:43.17 

0.0538 

4969, 4970, (575, 576) 

239 

Zw III 54 

iiizw54 

03:41:17.508 

+ 15:23:54.82 

0.0311 

4182 

23 

MKW 3S 

mkw3s 

15:21:51.708 

+07:42:24.65 

0.045 

900 

57 

MKW 8 

mkw8 

14:40:39.353 

+03:28:03.08 

0.027 

4942 

23 

PKS 0745-191 

pks0745- 

191 

07:47:31.265 

-19:17:41.62 

0.1028 

2427, 508, 6103 

56 

UGC 3957 

ugc3957 

07:40:58.133 

+55:25:38.25 

0.034 

8265 

8 

ZwCl 1215+0400 

zl215 

12:17:41.934 

+03:39:39.74 

0.075 

4184 

12 

ZwCl 1742+3306 

zl742 

17:44:14.447 

+32:59:29.02 

0.0757 

11708, 8267 

53 


a Cluster name. 

b Label used to denote cluster. 
c Cluster right ascention. 
d Cluster declination. 
e Cluster redshift. 

f The Chandra OBSID’s used in this work. OBSID’s in parantheses were used in imaging analysis, and excluded from spectral analysis. 
g The total exposure time for spectral analysis in kiloseconds. 
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acis_process_events. This reprocessing includes fil¬ 
tering to keep only event grades 0, 2, 3, 4 and 6. In 
addition, the VFAINT background cleaning method is 
applied to observations in VFAINT mode, using the 
check_vf_pha=yes option to acis_process_events. 

The center of the clusters is defined to be the centroid 
of event x- and y-coordinates, calculated using the fol¬ 
lowing iterative scheme. First, we select the observation 
with the longest exposure time, in the cases where we 
have multiple observations of one cluster. For ACIS-I 
pointing, we calculate the centroid using all four ACIS-I 
chips, while for ACIS-S pointings, we only use the chip 
with the most counts. We only include events with en¬ 
ergies between 0.3 and 7keV. For the first centroid com¬ 
putation iteration, Iteration 1, we calculate the medians 
and the standard deviations of the x- and y-coordinates 
of the events from the entirety of the selected chip(s) of 
the longest-exposure observation. For Iteration 2, we re¬ 
strict the median and standard deviation calculation to 
events within an ellipse with semi-major axes equal to 
twice the x and y standard deviations calculated in Iter¬ 
ation 1, i.e. 2 x (ct x i, cr.yi). For Iteration 3, the ellipse is 
shrunk to have semi-major axes equal to 1 x (cr X 2 , <J V 2 )- 
Similarly, Iteration 4’s filter ellipse has semi-major axes 
0.75 x (<J X 3 ,a y 3 ). For the final iteration, Iteration 5, we 
employ events from all observations in the centroid com¬ 
putation, instead of using the longest-exposure observa¬ 
tion alone. The Iteration 5 filter ellipse has semi-major 
axes equal to 1 x (a X 4 ,a y 4 ). 

Events from point sources are then identified and dis¬ 
carded. We use CIAO’s wavdetect tool, applied to an 
image of the merged event files from all observations. 
The input image to wavdetect includes only events with 
energies in the range 0.3-7keV and is binned in (2 x 2)- 
pixel bins. The detected point source regions are in¬ 
spected by eye to ensure that each region is large enough 
to include all events from its corresponding detected 
point source, and to add sources that were not detected 
by wavdetect. The latter tend to be point sources away 
from the telescope’s optical axis, where the point spread 
function is much larger than it is in the center. We also 
exclude any region of bright extended emission, which 
does not belong to the central cluster emission, such as 
that from infalling sub-clusters (e.g. the sub-cluster to 
the North of Abell 2163.) We expect that many clusters 
will contain emission from faint infalling sub-clusters, 
which cannot be resolved due to their low surface bright¬ 
ness. Therefore we do not attempt to discard all emission 
from identified sub-clusters, and only remove the bright 
peaks of such emission when present. 

Periods of high count rates resulting from flares are 
removed using the lc_clean() tool in Sherpa , which is 
CIAO’s tool for spectral analysis. We compute a light 
curve of all data counts, excluding the point sources de¬ 
tected above, and the central 300” to exclude the bulk 
of the cluster emission. Short flares are excluded when 
they are identified by lc_clean(), while longer flares are 
excluded manually by selecting events in the time range 
which is sufficiently removed from the the start or the 
end of the flare. Some observations in which one or many 
flares last for most of the exposure time are entirely ex¬ 
cluded. 

For each observation, we create a background dataset 
from the Blank Sky files, available as part of the Chan¬ 


dra calibration files. We choose the Blank Sky file for 
each ACIS chip based on the cluster dataset’s observa¬ 
tion date, its aim point and whether a CTI correction was 
applied to it. As mentioned above, we exclude all data 
taken prior to January 29, 2000. No Blank-Sky datasets 
are available in Period C, in VFAINT mode. We there¬ 
fore assign to these datasets, the Blank-Sky files from 
Period D in VFAINT mode. In addition, we exclude 
some ACIS-S data taken during Period C, because their 
corresponding Blank-Sky files are not available^. 

2.3. Data Analysis 

As is custom in clusters astrophysics, we define entropy 

as S = kT/nV ^, where k is Boltzmann’s constant, T the 
ICM temperature, and n e its electron density. As defined 
above, S is related to the thermodynamical entropy per 
particle, s, through s = (3k/2)hiS + sq, where so only 
depends on fundamental constants. We assume spherical 
symmetry and compute the entropy radial profiles from 
the density and temperature profiles, as described below. 

2.3.1. Density Profile, n e (r) 

Computing the electron density profile of the ICM is 
done in two steps. First, we use surface brightness (SB) 
measurements to constrain the shape of the density ra¬ 
dial profile. Then, we use spectral measurements to set 
the overall normalization of n e (r). To fit for the den¬ 
sity profile shape, we extract a SB profile based on pho¬ 
ton counts in the energy range 0.7-2keV. Radial bins are 
defined such that boundary radii are spaced logarithmi¬ 
cally, with a constant ratio of 1.25 between neighbor¬ 
ing radii, except when this spacing results in fewer than 
100 counts in the bin, in which case it is extended to 
the next radius. The innermost radius is defined as the 
projection of 1.2 arcseconds in the plane of the cluster, 
in order to include the largest number of counts, while 
avoiding any potential point source coinciding with the 
centroid of the ICM X-ray emission. From the source 
counts in each radial bin, we subtract the background 
contribution computed from the blank-sky datasets, and 
scaled by detector area and exposure time to match each 
observation and radial bin. The net number of counts is 
then normalized by the exposure map, to correct for the 
position dependence of Chandra's effective area. Finally 
we obtain a count SB by dividing by the solid angle of 
the extraction region. 

The emissivity of an X-ray plasma at cluster temper¬ 
atures is primarily in the form of bremsstrahlung and 
line emission. The contribution of each of those pro¬ 
cesses is proportional to n\, but depends differently on 
the plasma temperature. However, when emissivity is 
integrated over the 0.7-2keV energy range, the different 
temperature dependencies contrive to cancel each other, 
and the resulting emission in this energy band has a 
negligible temperature dependence. In addition, we as¬ 
sume the ICM to be optically thin. Therefore, the SB at 
any given point on the sky is simply the integral of all 
emission along the line of sight to that point. We per¬ 
form a maximum-likelihood fit to find the density profile 

6 The definition of the background periods used 
here is available in Maxim Markevitch’s note at 
http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 


2.3 Data Analysis 


5 


that best fits the counts SB profile, up to a normaliza¬ 
tion. We use an analytical form for n e (r ) that is flexible 
enough to allow us to phenomenologically fit the SB in 
all radial bins. Namely, we choose a simplified version o f 
Vikhlinin’s extended beta model (iVikhlinin et al.|[2006jh 


n 2 e {r) _ (£) _ 1 m 

n 2 0 (l + l|)3/3- a /2(l + (J_)3 )e /3 > ^ 

' c 

where no, r c , j3, a, r s and e are fit parameters. In each 
iteration of the fitting process, and for a given radial bin 
bound by radii r,; and 7 'i+i , we integrate « 2 (r) over the 
volume of the cylindrical shell defined by the above two 
radii and extending along the line of sight from —3Mpc to 
+3Mpc. The set of integrals from all radial bins is then 
compared to the corresponding SB values to determine 
the shape parameters in the right hand side of Equation 
□ 

This procedure provides the shape but not the nor¬ 
malization, no, of the density profiles. Computing no 
requires knowledge of the emission integral measure, 
El = f n e n p dV, where n p is the proton number density. 
This integral quantity, El, is simply proportional to the 
normalization of the APEC spectral model, whi ch we use 
to model ICM emission. The APEC model (ISnrith et al.l 
1200111 is fit to the spectra of multiple radial bins around 
the center of emission. The normalization of the spectral 
model of the i th radial bin, Ki, is related to Eli through 


0.82 10- 14 2ta 
^D\{1 + zf H ° 1 ’ 


( 2 ) 


where z is the cluster redshift and Da its angular di¬ 
ameter distance. We define Vj as the spatial integral of 
(n/no) 2 over the cylinder defined by the i th radial bin and 
bounded along the line of sight direction by l = ±20Mpc. 
We have assumed above that n e = 0.82n p , which is suit¬ 
able for typical ICM conditions. The best-fit no is deter¬ 
mined by minimizing 


2 _ (JC ^ CpVjnl) 

Wo / j SK 2 ^ 

i 1 

with respect to 71 q, where SKi is the uncertainty on K, 
and C D = 0.82 10~ 14 /ATrD 2 A (l + z) 2 . 

2.3.2. Temperature Profile, kT(r) 

To compute the temperature profiles, we again con¬ 
struct radial bins, and fit their spectra, resulting in a pro¬ 
jected temperature profile. We then deproject the above 
temperature profile, using similar methods to those used 
to deproject the density profile. We choose the size of 
the radial bins in order to include the lowest number of 
counts necessary for a temperature determination, with 
10% uncertainty. This count number, N nec , which is 
a function of temperature and background count frac¬ 
tion, is estimated by simulating spectra of different total 
counts, temperatures and background fractions, f bs , to 
find the n ecessary coun ts for a 10% uncertainty on tem¬ 
perature (lElkholvll2012D . We find that N nec can approx¬ 
imately be fit by 


/ hj 1 \ 1 - 7 

N nec (kT, f bg ) = 500 x 10 1 - 976 *. x l -—j . (4) 

We extract spectra of at least N nec counts and fit their 
spectra to an absorbed APEC model in CIAO’s tool, 
Sherpa. We create a spectrum from each observation that 
partially or wholly covers the annulus corresponding to a 
radial bin, and simultaneously fit these spectra. We link 
the temperatures and metal abundances of these spectra, 
across the multiple observations, during the fit. Normal¬ 
izations are only linked for spectra that cover more than 
95% of the solid angle of the annular region. Metallicity 
is left as a free parameter, while the hydrogen column 
de nsity, nn, is fix e d to t he values from the LAB dataset 
of iKalberla et al.l (12005!) . We also incl ude i n our o nline 
datase t the results of analysis with iDickev fe Lockmanl 
(H99l nn measurements. We find that leaving nn as 
a free parameter returns unreasonable best-fit values on 
both nn and temperature. However, for Abell 478 , which 
is reported to have varying nn bv lYikhlinin ct al.l (12005 1. 
we allow nn to be free, with a minimum equal to the 21- 
cm-measured value. 

The background spectrum is extracted from the Chan¬ 
dra blank sky datasets and from source free regions 
around cluster observations. The background is modeled 
as particle background, plus X-ray background. The X- 
ray background is modeled as an absorbed 0.2keV APEC 
model, with nn = 2.09 x 10 22 cm -2 and an absorbed 
power-law component for the Cosmic X-ray Background, 
with index set to -1.4. These components are convolved 
with the instrument response for each CCD chip. The 
particle background is modeled as a series of Gaussian, 
exponential, and a power-law functions, to phenomeno¬ 
logically fit the remaining components of the blank sky 
datasets. It is not convolved with the instrument re¬ 
sponse. The overall background spectral model varies 
from epoch to epoch, and also depends on the CCD chip 
used. 

We assume that the shape of the instrumental back¬ 
ground component of the cluster observation is the same 
as the best-fit model from the corresponding blank sky 
dataset and compute its normalization as suggested in 
Maxim Markevitch’s cookbook for treating the back¬ 
ground datEQ The overall normalization of the instru¬ 
mental background is computed by scaling the back¬ 
ground normalization according to the ratio of counts 
in the 9.5-12 keV energy range, in the cluster dataset 
relative to the blank sky dataset. In this manner we at¬ 
tempt to capture any possible change in the background 
normalization between different epochs. We note that 
datasets in our sample with OBSID between 7686 and 
7701 are missing high-energy counts. Their instrumen¬ 
tal background normalization is thus scaled simply by 
exposure time and solid angle. 

Having used the blank sky data to constrain the instru¬ 
mental background components, we proceed to fit the X- 
ray background from the in-field spectra. The latter are 
modeled with the same model described above plus an 
additional APEC component to account for residual clus¬ 
ter emission. The APEC model’s temperature is fixed at 

7 http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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2 DATA SAMPLE AND ANALYSIS 


the temperature measured outside a projected radius of 
150kpc, using an initial simple fit. 

The in-field spectra are obtained from annuli centered 
around the cluster center, and covering regions that are 
visually identified to contain mostly background X-ray 
counts. For Abell 119 and Abell 3571, the cluster emis¬ 
sion covers most of the field of view (FOV.) We thus rely 
on the blank sky background data to model both instru¬ 
mental and X-ray background components, for these two 
clusters. 

We use CSTAT as our fit statistic, as it is more suitable 
for energy bins with low counts, where the more com¬ 
monly used x 2 statistic introduces bias. After obtaining 
the best-fit temperature, we compute its uncertainty us¬ 
ing Sherpa’s proj () function, which varies temperature 
along a grid and searches for the best-fit at each temper¬ 
ature by varying the other thawed parameters. 

The spectral fitting described above, returns a best-fit 
projected temperature for a given radial bin: Since the 
ICM is thought to be optically thin, the emission at one 
point on the sky is the sum of all emission from the line of 
sight behind that point. Thus, to compute the true three- 
dimensional temperature profile, we assume a flexible an¬ 
alytic form for kT(r), vary its parameters repeatedly, 
projecting it along the line of sight in each iteration un¬ 
til the best match is found with the measured projected 
temperature radial profile. This fitting process is again 
run using Sherpa. Th e projection is compute d according 
to the prescription in iMazzotta ct akl (12004 1. who show 
that to recover a single-temperature fit from a mixture of 
many temperature components, one should average these 
temperatures with a weighting proportional to n 2 VT a , 
where V is the volume of the region of emission. We 
choose a = —0.75, as suggested by the range of values 
found bv IMazzotta et alT(l2004 1 for spectra of different 
metallicities. Th e three-dimensional te mperature profile 
is modeled as in IVikhlinin et akl (120061) : 


kT{r) = kT 0 


( r/r t ) a 
[1 + (r/r t ) b } c/b 


x + T m i n /TQ 
x + 1 


( 5 ) 


where x = (r/r coo ;) a =° oi and where T 0 , r t , a, b , c, T min , 
vcool and a coo i are fit parameters. The number of free 
parameters depends on the number of available temper¬ 
ature measurements. 

Computing kT(r) and n e {r) gives us the necessary 
quantities to measure i?soo, the gas mass within i? 5 oo 5 
which we call M gas , and the total gravitatio nal mass 
withi n the same radius, As described in lElkhoTvl 

(Hm, we do so using an iterative scheme, since the 3 
quantities are related . We use the M — Yx relation of 
iKravtsov et al.l (120061 ) to relate M 500 to our measurables, 
as Yx = kT x M gas . 


2.3.3. Uncertainties 

The uncertainties on the entropy profile of each clus¬ 
ter are estimated by generating a set of S(r) models, 
which are allowed by the data and their uncertainties, 
as descibed here. For temperature radial profiles, using 
the temperature measurements in each radial bin, and 
their error estimate, we randomly generate new “fake” 
datasets, and fit them one at a time. To generate a fake 
temperature measurement for each radial bin, we draw 


its value from a random distribution designed to capture 
the asymmetric uncertainties obtained on the bin’s best- 
fit temperature. This probability distribution is a piece- 
wise function of 2 Gaussian distributions on either side 
of the best-fit temperature, with the standard deviations 
equal to the measured 1 -a upper and lower uncertainties. 
The latter are not in general equal to each other. Once 
a complete radial temperature profile is generated over 
the entire available radial range, we fit it with the same 
model in Equations 0 and repeat this analysis for 400 
iterations. 

The same analysis is repeated for the density profiles, 
where surface brightness measurements are similarly per¬ 
turbed for 300 iterations according to their uncertainties. 
The uncertainties in this case are assumed to be symmet¬ 
ric, and the fake surface brightness measurements are 
drawn from a Gaussian distribution. 

To compute the uncertainty on the entropy profile, 
we compute a set of entropy profiles from pairing dif¬ 
ferent temperature and density profiles, from the above 
generated sets. We iterate through all 400 tempera¬ 
ture profiles. For each temperature profile, we iterate 
through 10 density profiles computing a temperature pro¬ 
file S(r) = kT(r)/n e (r ) 2 / 3 , in each iteration. We ensure 
to choose different density profiles, from one tempera¬ 
ture profile to the next, until all 300 profiles are used, at 
which point we start from the beginning of the density 
profiles list. The result is an ensemble of 4000 entropy 
profiles, which we use to find the distribution of entropy 
values at any given radius. 


2.3.4. Metallicity Profile Calculation Method 

To compute the metallicity profiles, we take a simi¬ 
lar approach to that used to make the temperature pro¬ 
files. First, using simulated spectra, we estimate the 
minimum necessary counts, N^ ec (kT 1 Z, fbg), to obtain 
a 20% uncertainty on the best-fit metallicity. In this 
case, N% ec (kT, Z, fb g ) does not have a simple analytical 
form as its counterpart for temperature measurement, 
but is rather estimated from a weighted average of N% ec 
values estimated for the fcT, Z and fb g values that were 
simulated (lElkholvll2012D . Then, using the derived N^ ec , 
we extract spectra in radial bins using the same bins 
used fo the kT(r ) profile calculation, and joining them 
whenever more counts are needed for a 20%-uncertainty 
temperature estimate. We take the maximum radius of 
extraction to be i? 5 oo- 

For spectral fitting, we again model both source emis¬ 
sion from the ICM and background. For the background 
spectra, we use the same best-fit parameters found in 
the kT(r) analysis, above. The background normaliza¬ 
tion is computed using the same method as in the kT(r) 
analysis, described in Section f2.3.21 Cluster emission is 
modeled using a 1-temperature model, and using a 2- 
temperature model, where the cooler component’s tem¬ 
perature is set to one half of the value of the hotter com¬ 
ponent’s temperature. 

We compute the uncertainty on the metallicity in each 
bin using Sherpa ’’s projO function. From the obtained 
metallicity profile, we characterize the metallicity of each 
cluster by two global quantities, Z m a and Z ln . We define 


Ymid, — 


So.l5<n<0.3_R 5 oo ZjMgas^i 
So. 15<n<0.3fl 5 oo M g as,i 


( 6 ) 
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which is the gas-mass-weighted metallicity over all shells 
in the range 0.15 < r < 0.3i?5oo- Here, Z,j and M gaSti 
are, respectively, the metallicity^ and the gas mass in the 
i th radial bin. In other words, Z m id traces the total iron 
mass, in the region 0.15i?5oo < r < O. 3 .R 500 , ac¬ 

cording to 

= A Fe Z mid M™ d , (7) 


where p and p e are the mean number of nucleons per par¬ 
ticle and per electron, respectively, m p is the proton mass 
and p c {z) is the Universe’s critical density at the redshift 
of observation, z. The characteristic entropy, 6*500 at an 
overdensity S = 500 is simply the entropy obtained using 
the characteristic temperature at S = 500, which is the 
equivalent of the virial temperature but defined for R 500 
instead of the virial radius, 


where Ap e = 0.0019 is the solar abundance of iron by 
mass, accordin g to the photospheric measurements in 
lAnders fe Grevessel (119891) assumed for our spectral anal¬ 
ysis, and M™ d is the gas mass contained in the same 
region. Similarly, 


Er-j<0-15^500 ZjMg aSt i 
Sri<0.15fl 5 oo Mgas,i 


( 8 ) 


Hereafter, Z m i d will be used as a measure of the metal¬ 
licity of the bulk of a cluster, while Zi n will be used as a 
measure of the core metallicity in clusters. 


3. RESULTS 

We make our data available on the FTP sitc@ 


ftp://space.mit.edu/pub/tamer/ebc2015/ . 

Appendix O contains the plot of the entropy profile 
for each cluster, while the individual metallicity radial 
profiles are shown in Appendix [Bl 

We first note that 4 clusters of our sample deviate 
considerably from spherical symmetry simply based on 
their surface brightness image. These clusters, Abell 754, 
Abell 2256, Abell 3376 and Abell 3667, are known to be 
undergoing merging events. We exclude the disturbed 
clusters from our analysis, except when noted. 

Table ICTl in Appendix IH1 shows the measured values of 
kTx , R 500 , A /500 an d M gas for each cluster. 

3.1. Entropy Profiles 

Our best-fit results for density and temperature radial 
profiles are shown in Tables [DTI and ID21 respectively, in 
Appendix [D] 

We overplot the computed entropy profiles for all of 
our sample’s clusters in Figure [Q On the left panel, 
we plot entropy as a function of radius, which we nor¬ 
malize with respect to i?soo- On the right panel, we 
plot entropy as a function of enclosed gas mass frac¬ 
tion, F g = M g /(fi,M 500 ), where M g is the interior gas 
mass, and fb the universal baryon fraction with respect 
to all matter, i.e. fb = flb/fl m - We use F g to plot 
entropy profiles because this is the variable used in a 
Lagrangian des cription of the entropy distribution in 
c lusters (see e.g. iTozzi fe Normanll200ll fVoit et alJl2003t 

iNath fe Maiumdail 1 201 Ilk The en tropy in both p anels 
is normalized with respect to 6500 (IVoit et alJl2003 ) , the 
characteristic entropy of the cluster at R 500 : 

_ GM 500 pm p 

>3500 =-577 > (,yj 

2i?500 [500 fbpc{z)/PerripY' 

8 The file named “README” within this FTP site details the 
content of the data. 


kTsoo 


GM 500 pm p 

2i?500 


( 10 ) 


and using the average electron density inside R 500 , 


n e = 500 fbp c (z)/p e m p . (11) 

It represents the entropy scale set by gravity in the self¬ 
similar picture. 

The dark blue line in Figure [T] represents a model of the 
entropy profile of a cluster generated from gravitational 
collapse alone, which w as calculat ed with hydrodynami- 
cal AMR simulations in lVoit et all (j2 005l ). Voit’s entropy 
profile is approximated analytically as a power law, 

/ \ 1.24 

Sv(r) = 1.536500 () , (12) 

V ^500 / 

and is valid for radii larger than ap proximately 0 . 26 * 500 . 
We em ploy the conversion used in iPointecouteau et al.l 
(|2005l f to express Equation [12] in terms of 6500 and R 500 , 
as opposed t o the measuremen ts at an overdensity of 200 , 
presented in lVoit et al.l (120051 b 

Turquoise curves in Figure [l] represent cool core clus¬ 
ters (CC), while red curves are for non-cool core clus¬ 
ters (NCC). We use th e surface br igh tnes s concentra¬ 
tion, csBi introduced in lSantos et al.l (I2008D to quantify 
the cool core state of a cluster. The parameter csb is 
defined as the ratio of the surface brightness within 40 
kpc to that within 400kpc of the pea k of the emissior0. 
The values of Csb are shown in Table lEll We define CC 
clusters as clusters with csb > 11, while NCC cluster 
have csb < H- 

The first observation to make is that for most of the 
studied radial range, all entropy profiles lie above Voit’s 
2005 gravitationally induced entropy m odel. T his re¬ 
sult has been known in the literature (e.g. lPonman et all 
Il999t iPratt et aDI l2010ll . and the additional entropy in 
the observations has been attributed to non-gravitational 
processes, such as winds and AGN heating. Second, 
some of the entropy profiles in Figure [0 show a de¬ 
crease or a flattening starting at a radius between 0.3 
and 0 . 8 i? 5 oo- This is due to a decrease in measured tem¬ 
perature towards the outskirts of many clusters, which is 
not matched by a steep enough decrease of measured den¬ 
sity with radius. Such a configuration of large amounts 
of lower-entropy gas at larger radii is not physically sta¬ 
ble. There are two sources of systematic error that can 
be contributing here to give erroneous temperature and 
density measurements. First, systematics in our estimate 
of the level of X-ray background will translate into an 
error in the estimate of the cluster surface brightness at 

9 For calculating csb , we use the emission peak as the clus¬ 
ter center, unlike in previous analysis where the emission is used 
instead. 






































3 RESULTS 




FIG. 1: Plot of all computed entropy profiles. Left: scaled entropy, S/S 500 , as a function of scaled radius, r/Rsoo■ 
Right: scaled entropy as a function of enclosed gas mass fraction, F g . Turquoise curves are for CC clusters, while red 
curves are for NCC clusters. The dark blue curve is the power law describing the Voit 2005 entropy profile found in 
hydrodynamical simulations. Dashed curves are for the 4 morphologically disturbed clusters Abell 754, Abell 2256, 
Abell 3376 and Abell 3667. 


these large radii, introducing a bias to the inferred den¬ 
sity, and also biasing outer temp erature estimates. Sec - 
ond, the deprojection method of lMazzotta et al.l (120041) . 
which we employ to deproject the measured 2D temper¬ 
ature profile to a 3D kT(r), is known to be less accu¬ 
rate when there is significant contribution to the emis¬ 
sion from spec tral components with t emperatures sm aller 
than ~ 3keV dMazzotta et al.ll2004t IVikhlininl 120061) . 

In addition, we show in Figure ET1 in Appendix [Al the 
individual entropy profiles we compute. The estimate 
of the 1-sigma range of entropy at each radius is repre¬ 
sented by the turquoise shaded region. The uncertainty 
in the measured temperatures, which translates into an 
uncertainty in the parameters of the temperature radial 
profile, is the main contributor to the uncertainty in the 
entropy profile. By comparison, the contribution of the 
density uncertainty to the entropy uncertainty is much 
smaller. 

The light gray error bars, in Figure EH represent 
th e entropy profiles meas ured in the ACCEPT study 
by iCavagnolo et al.l (12009D . Our entropy profiles agree, 
in general, with the ACCEPT entropy profiles, where 
they overlap. However, we extend our entropy profiles to 
larger radii, where we model both density and tempera¬ 
ture. 

3.2. Metallicity Profiles and Global Measurements 
3.2.1. Profiles 

We overplot all obtained metallicity profiles in Fig¬ 
ure [2] The dispersion in the values of observed cluster 
metallicities decreases for radii larger than 0.1 — O. 2 .R 500 , 
despite the larger uncertainty associated with measured 
metallicities at these high radii. One metallicity mea¬ 
surement seems to be exceptionally larger, at large ra¬ 
dius, as seen in Figure [2j This is the last metallicity 
measurement for Abell 2204, where modeling the back¬ 


ground is likely suffering from systematics, despite an 
acceptable fit statistic. The outermost spectra for Abell 
2204 are found to have excess low-energy counts, which 
are mostly fit by our galactic X-ray background compo¬ 
nent. 



r/, H500 

FIG. 2: Superimposed metallicity profiles of clus¬ 
ters. Red curves are for NCC clusters, while 
turquoise curves are for CC clusters. 

We also present in Appendix [5] plots of our measured 
metallicity profiles, individually, with their esti mated un¬ 
certainty. The two vertical lines in Figure IB 11 are drawn 
at radii O.I 5 .R 500 and I? 5 oo, for those clusters where we 
can measure M 500 and R 500 using our iterative method. 

3.2.2. Core Metallicity, Zi n 
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The measured values of Zi n are shown in Table ED in 
Appendix [F] We plot Z in as a function of cluster mass, 
and as a function of the global temperature measure¬ 
ment, kTx , in Figured Figure [3] suggests that low-mass 
clusters ( M 500 < 3.5 x 10 14 Mq) exhibit more dispersion 
in Zi n than high-mass clusters do. 

To quantify this latter observation, we show in Table 
[2] the error-weighted mean, and the standard deviation 
of the values of Z in . for clusters grouped by mass. The 
results for the 1-temperature fit are labeled “IT”, and 
those for the 2-temperature fit “2T.” The standard de¬ 
viation of Zi n for clusters with M$oo > 3.5 x 1O 14 M 0 
is 2.4 to 2.9 times smaller than that for clusters with 
-^500 < 3.5 x 10 14 Mq, for the IT and the 2T models, re¬ 
spectively. The value of \ 2 for each subset shows whether 
the dispersion for low-mass clusters is solely the result of 
measurement uncertainties. Here, 

2 _ (Zin,i - (Zj n )) , . 

x 2^ sz 2 . ’ [ ’ 

i in,i 

where the sum is over the sample of clusters denoted by 
i, Zin,i is the i th cluster’s central metallicity, <5 Zi n ,i is its 
uncertainty and ( Zi n } the sample mean. 

For low-mass clusters, we find that y 2 is more than 30 
times the number of degrees of freedom, for both IT and 
2T fits. This confirms that the dispersion seen in low- 
mass clusters — in the range (0.12 — 0.15 )Zq — is not 
driven by measurement uncertainties. On the contrary, 
the dispersion of Zi n in high-mass clusters has a more 
significant contribution from measurement uncertainties, 
despite being much smaller at ~ 0.05Zq. 


TABLE 2: Statistics for Zi n . 



Model 

Mean 

Std. Dev. 

X 2 /dof 



(Zq) 

(Zq) 

M 50 o < 3.5 x 10 14 M 0 

IT 

2T 

0.30 

0.25 

0.12 

0.15 

33. (490./15) 
31. (470./15) 

M 500 > 3.5 x 1O 14 M 0 

IT 

2T 

0.36 

0.35 

0.050 

0.052 

2.6 (64,/25) 
1.2 (29./25) 


Note. — Error-weighted mean, standard deviation of Zi n and 
X 2 /dof with respect to calculated mean. Results shown for clusters 
with M 500 smaller than and larger than 3.5 x 10 14 Mq, and for fits 
using the IT and the 2T spectral models 


The above results are unchanged when we include 
the Zi n measurements from the high-mass, asymmetrical 
clusters Abell 754, Abell 2256 and Abell 3667, which have 

Zin = 0.369±°;°35 , 0.389t°;°6i and 0.345 re¬ 
spectively. These clusters were excluded from the above 
analysis based on their asymmetric morphologies. The 
IT dispersion of Zi n , over the high-mass subset decreases 
to 0.047Zq, when the above 3 clusters are included, while 
the 2T dispersion does not change significantly. 

To emphasize the difference between low- and high- 
mass clusters metallicity dispersions, we show in Figured 
the superimposed metallicity radial profiles of all clusters 
in our sample, where we differentiate between the two 
subsets by color. Green data points represent low-mass 
clusters, while orange data points represent larger clus¬ 
ters. We can see in this plot that large clusters’ metal¬ 


licity profiles are less scattered than low-mass clusters’. 
This translates in the different observed dispersions of 
Z^ seen in Figure [3l 

Comparing CC and NCC clusters, we find that for CC 
clusters Z in = 0.30 ± O.IOZq, while for NCC clusters 
Z in = 0.37 ± 0.080Z q . We find that Z in is larger for 
NCC clusters compared to CC clusters in this sample, 
despite the difference being within the measured disper¬ 
sions of both quantities. This contrasts to the metal 
excess measured in t he cen ters of CC cluster in e.g. 
IDc Grandi fe Molcndil (120Q1H . 

3.2.3. Outer Metallicity, Zmid 

The measured values of Z m id are shown in Table iFlT in 
Appendix IF1 In Figured we show a plot of Z m id , which 
measures the average metallicity outside the core, as a 
function of the total mass, M 500 , and as a function of 
kTx ■ It Tabic [3] we show the statistics for the distribu¬ 
tion of Z m id values. In the case of outer metallicity, we no 
longer detect a clear difference in the dispersions of high- 
and low-mass clusters. However, when we compare Z m id 
to Z ^, we find that for each of the low- and high-mass 
cluster samples, Z m id values are smaller than Z ln values. 
For high-mass clusters, for example, the mean and stan¬ 
dard deviation for Z m id are (0.27 ± 0.073 )Zq, while for 
Z in they are (0.36 ± 0.050 )Zq. This points to a decrease 
in the iron mass fraction as we move from the core region, 
r < 0.15i?5oo, to the outer region 0.15 < r < 0 . 3 i? 5 oo- 
This decrease_is however within the measured dispersions 
of Z m id and Zi n and is also found for low-mass clusters. 


TABLE 3: Statistics for Zmid- 



Model 

Mean 

Std. Dev. 

X 2 /dof 



(Zq) 

(Zq) 

Mbqo < 3.5 x 10 14 M o 

IT 

2T 

0.22 

0.20 

0.079 

0.085 

10. (150./15) 
9.1 (140./15) 

M 500 > 3.5 x 10 14 M q 

IT 

2T 

0.27 

0.27 

0.073 

0.074 

3.2 (79./25) 

2.3 (56./25) 


Note. — Error-weighted mean, standard deviation of Z m id and 
X 2 /dof with respect to calculated mean. Results shown for clusters 
with M 500 smaller than and larger than 3.5 X 10 14 Mq, and for fits 
using the IT and the 2T spectral models 


3.3. Metallicity-Entropy Relation 

One direct approach to look for a relation between 
pre-enrichment and pre-heating is to look for a corre¬ 
lation between the ICM non-gravitational entropy and 
the ICM bulk metallicity measured outside the central 
region of the cluster, Z m id • In addition, we also con¬ 
sider the relation between non-gravitational entropy and 
Zi n . We use the ratio of measured entropy to the ex¬ 
pected gravitational entropy, S grav , to probe the amount 
non-gravitational entropy. We define the scaled entropy 
x s = S/S grav , which we use as a measure for any non- 
gravitational entropy, and discuss our assumptions on 
ffgrav? below. 

The metallicity-entropy relation is studied with en¬ 
tropy measured at several locations in the clusters. 
First, we consider entropy measurements at fixed R 500 - 
scaled radii. This is justified because the gravita- 
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3 RESULTS 



FIG. 3: Core metallicity as a function of cluster mass or temperature. Left: Zi n vs. Msoo- Right: Zi n vs. kTx . Red 
points are for NCC clusters, while turquoise points are for CC clusters. 
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FIG. 4: Superimposed metallicity profiles of clus¬ 
ters. Orange data points are for clusters with mass 
A /500 > 3.5 x 1O 14 M 0 , while green data points cor¬ 
respond to clusters with JV /500 < 3.5 x 10 14 Mq 

tional entropy model of IVoit et ahl (120051 ) scales self- 
similarly, and is given in terms of a profile which is a 
function of r/R^o- In this case, we define S grav as 
the expected gravitational entropy from Voit’s model, 
Sgrav (r) = S v (r). Second, we take a Lagrangian ap¬ 
proach and study entropy at a fixed interior gas mass 
fraction, F g = M„/(fhMso 0 ), where M„ is the interior 
gas mass (see e.g. iTozzi fc Normanll2001l : IVoit et alJ2003t 

[Nath fe Maiumdarll2011lh ~ This can be useful because 

buoyancy tends to order the ICM such that low-entropy 
gas finds its way to the bottom of the cluster potential, 
while high-entropy gas rises to large radii. In this latter 
case, we simply use S grav = Ssoo to scale the entropy, 
to avoid using a specific model of entropy dependence 
on F g , while still capturing the 5500 scaling expected in 



M 500 (M s ) kT x (keV) 


FIG. 5: Outer metallicity as a function of cluster 
mass or temperature. Left: Zmid vs. Msoo- Right: 

Zmid. vs. kTx . Red points are for NCC clusters, 
while turquoise points are for CC clusters. 

self-similar galaxy clusters. 

In the analysis below not all clusters are included for 
each measurement. The first filter we apply is to exclude 
4 clusters, which are visually judged to greatly deviate 
from spherical symmetry. These are Abell 754, Abell 
2256, Abell 3376 and Abell 3667. In addition, in the suc¬ 
cessive measurements at different radii, below, we only 
include a cluster at a certain radius if the size of the FOV 
is larger than the radius of interest. 

We start by looking at the Z m id — x s relation at con¬ 
stant scaled radius. As described above, in this case, 
x s = x 3 (r) = S (r)/ Sg ra v (r) = S(r)/S v {r). We mea¬ 
sure entropy at r = 0.2R 500 , 0.3i? 50 o, 0.5i? 5O o, O.8i? 50 o 
and li? 5 oo- Our measurements are shown in Figure [6] 
Similarly, we consider the same relation at fixed F g . 
We choose values of F g corresponding to the sample 
average across all clusters, at r = 0.2, 0.3, 0.5, 0.8 
and li? 5 oo- Table [4] shows the correspondence between 
scaled radius and the sample average gas mass fraction. 
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Our metallicity-entropy measurements at constant F g are 
shown in Figure [7] 

TABLE 4: Cluster sample average of F g at various 
radii. _ 


r/Rsoo 

0.2 

0.3 

0.5 

0.8 

1.0 

F, 

0.050 

0.10 

0.24 

0.47 

0.63 


As can be seen in Figures [5] and [3 there is no visible 
correlation between our estimate of non-gravitational en¬ 
tropy, and_the bulk metal content of a cluster, as esti¬ 
mated by Z mid . 

As for metallicity measured in the core, Zi n , we expect 
low-radius metallicity measurements to probe processes 
that occur after the collapse of the yluster. We repeat 
the analysis performed above, with Zi n instead of Z m i d . 
Figure [5] shows plots of the measured inner metallicity, 
Zi n , against the ratio of measured entropy to Sy, at the 
above-mentioned scaled radii. Again, there is only weak 
indication of a correlation between inner metallicity and 
x s , at smaller radii. We perform statistical analysis using 
a bootstrap resampling method to calculate the signifi¬ 
cance of the correlation between the various metallicity 
and entropy measures. The lowest obtained p-values are 
of 1.3% and 2.0% for the CC-only samples at r = 0.3 and 
O. 2 .R 500 , respectively. 

4. DISCUSSION 

4.1. Inner Metallicity Scatter Systematics 

As was presented in Section f3.2.21 the core iron mass 
fraction, Zi n , over our cluster sample has a different dis¬ 
tribution for low-mass clusters than for high-mass clus¬ 
ters. Measurements of Z in in large clusters ( M > 3.5 x 
10 14 Mq) are narrowly distributed around their mean of 
0.36 Zq, with a standard deviation of only oz = 0.050 Zq. 
On the other hand, Zi n for low-mass clusters has a stan¬ 
dard deviation of oz = 0.12 Zq, around a slightly lower 
mean value for the sample. See Table [2] for details. 
The uncertainties on the individual Z.- m measurements 
are too small to explain the dispersion in low-mass clus¬ 
ters, since x 2 /d.o.f = 33. This means that the observed 
scatter is intrinsic to the data, and not a result of mea¬ 
surement uncertainties. On the other hand, we calculate 
X 2 /d.o.f = 2.6 for Z in in large clusters, indicating that 
measurement uncertainties contribute relatively more to 
the scatter, which nonetheless has a much lower value of 
only oz = 0.050 Zq. 

We attempt here to understand the difference between 
the distribution of low- and high-mass clusters’ inner 
metallicity values. First we check whether the observed 
effect is due to systematics, and then we present sev¬ 
eral physical explanations of the measurements, in the 
following sections. 

The first systematic effect to be tested concerns the 
inclusion of X-ray photons from the region r > 0.15i?soo 
in the computation of Z; m . This occurs because spectral 
regions for metallicity measurement are_defined before 
i? 50 o is computed, while we desire for Zi n to measure 
the metallicity within 0.15i?5oo- Once R 500 is calculated, 
to compute Zi n , we include all radial bins which overlap 
the disc r < 0.15i?soo. The last such bin will, in general, 


extend beyond r = 0.15i?5oo- We address this by esti¬ 
mating the fraction of counts originating from the region 
r > 0.15i?5oo, which are used to compute Zi n . We de¬ 
note this count fraction by f ou t, and show a histogram 
of its distribution in our sample in Figure [TUI 

Since metallicity profiles generally decrease with ra¬ 
dius, including emission from large radii for a given clus¬ 
ter might bias its Z, in measurement, compared to the 
rest of the sample. We would like to test the magnitude 
of this effect. We thus repeat the measurement of Zi n 
dispersion with a sample which excludes clusters with 
significant contribution from r > 0.15i?soo. We choose 
the cutoff value to be f out = 0.25 to capture the peak of 
the distribution of clusters with f out < 0.25 seen in Fig¬ 
ure [10] We calculate the mean, dispersion and x/d.o.f 
statistics as was done with the complete sample. The 
results are displayed in Table [5] We find that the differ¬ 
ence between the high- and low-mass cluster samples still 
remains for dispersion and X/d.o.f. We thus conclude 
that the spectral bins’ sizes do not have a large effect on 
this discrepancy in measured metallicity dispersions. 


TABLE 5: Statistics for Zt n excluding clusters 
with f ou t > 0.25. 



Model 

Mean 

Std. Dev. 

X 2 /dof 



{Zq) 

{Zq) 

M 50 o < 3.5 x 10 14 Mq 

IT 

2T 

0.30 

0.26 

0.15 

0.18 

59. (470./8) 
57. (460./8) 

M s oo > 3.5 x 10 14 M q 

IT 

2T 

0.36 

0.35 

0.038 

0.043 

3.7 (49./13) 
1.5 (19./13) 


Note. — Similar to Tablc[2] using a cluster sample that excludes 
fout > 0.25 clusters. 


The second systematic effect we test is the effect of 
the number of radial bins used to measure Zi n in the 
obtained value. As can be seen in Equation [8j Zi n is 
a weighted sum of single metallicity measurements. For 
clusters without enough photons to create multiple ra¬ 
dial bins within 0.15i?5ooj the measurements will give 
less precise estimates of Z in , on average. We find^ how¬ 
ever, that there is no significant dependence of Zi n on 
the number of bins used to estimate it. In addition, we 
also find that the number of bins covering r < O.I 5 .R 500 
does not depend on M 500 . 

We also repeat the analysis using a constant physical 
radius aperture of 150kpc to compute the inner iron mass 
fraction. We find that the different levels of dispersion 
remain unchanged, even with the physical radius aper¬ 
ture. 

4.2. Metallicity Scatter as a Reflection of Structure 
Formation 

One possible physical explanation to this observed dif¬ 
ference between low-mass and high-mass clusters, is that 
the metal content in clusters is driven by the merger his¬ 
tory of clusters. In the hierarchical model of structure 
formation, low-mass clusters, groups and galaxies merge 
to form the larger-mass clusters. Thus, if the metal con¬ 
tent is non-uniform across all these progenitors, as they 
merge with each other, the resulting metallicity is an av¬ 
erage of the initial progenitor metallicities. In a simple 
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FIG. 6: Plot of gas mass-weighted metallicity Zmid , vs. the scaled entropy, S/Sv■ Panels from left to right represent 
the cases where entropy is measured at r = 0.27?soo, r = 0.3-Rsoo, r = 0.5f?soo, r = O. 8 -R 500 and r = Rsoo- Only 
clusters with Chandra coverage at each of these radii are represented in the corresponding panel. Turquoise data points 
are for CC cluster, while red points represent NCC clusters. Abell 400 is the outlying cluster with large metallicity at 
Zmid ~ 0.6Zq. 



FIG. 7: Plot of gas mass-weighted metallicity Z rnl( i, vs. the scaled entropy, S/Sv■ Panels from left to right represent 
the cases where entropy is measured at F g = 0.050, F g = 0.10, F g = 0.24, F g = 0.47 and F g = 0.63. Only clusters 
with Chandra coverage at each of these radii are represented in the corresponding panel. Turquoise data points are 
for CC cluster, while red points represent NCC clusters. Abell 400 is the outlying cluster with large metallicity at 
Zmid ~ 0.6 Zq. 


model, if the metallicities of these progenitor structures 
are distributed around a mean universal value, Zq, then 
the sum metal content in a cluster formed by the merger 
of all these components should approach Zq, as the num¬ 
ber of components increases. The iron mass fraction, Zi n , 
of a large cluster will thus be an average of the metallic- 
ities of its smaller progenitors. 

Under this hypothesis, the decrease of the dispersion 
of Zi n as we go from low-mass clusters to high-mass clus¬ 
ters simply results from the mixing of low-mass clusters’ 
gas, after they merge to make larger clusters. The mixing 
then results in averaged, less dispersed metallicity values 
in the merged clusters, compared to the initial progeni¬ 
tors’ metallicities. It must be noted, however, that a high 
mass cluster from our sample, say of M 500 ~ 1O 15 M 0 , 
will not be exclusively formed by the merger of 1O 14 M 0 - 
sized clusters, i.e from the low-mass extreme of our sam¬ 
ple. A 1O 15 M 0 cluster will undergo numerous mergers 
involving galaxy- and group-sized haloe s, as well as a 
smoo th and continuous accretion (see e.g. lFakhouri et all 
I20l0lh This does not contradict the above hypothesis, as 
the more numerous the components making up a cluster 
are, the closer its metallicity approaches the universal 


average value. 

If this is the correct explanation for the observed larger 
dispersion of Zi n in low-mass clusters, then_it should also 
be reflected in the outer radii metallicity, Z m idi which is 
measured in the range 0.15 < r < 0 . 3 i? 5 oo- Our analy¬ 
sis however does not detect the same signal in the outer 
regions, as we do_ for inner metal content. The sam¬ 
ple dispersion in Z m id is roughly 0.079 Zq (O.O73Z 0 ) for 
the low-mass (high-mass) sample, and the contribution 
of measurement uncertainties to that dispersion is esti¬ 
mated to be around O.O4i ? 0 (0.05 Zq) for the low-mass 
(high-mass) sample. 

The decrease of the dispersion of inner metallicity from 
small to large clusters requires that, as clusters merge, 
metals from the progenitor clusters are able to efficiently 
find their way to the center of the cluster, while avoiding 
mixing with gas in the outer regions of clu sters. Inter¬ 
estingly, the hydrodynamical simulations of lCoral (2006) 
credit the infall of cold metal-rich clumps from large radii 
for the metal enrichment of the cluster central regions. 
A similar p roce s s is a lso one of the mechanisms evoked 
in iMillion et all ( 2011 1) to explain the existence of a peak 
in the observed radial profiles of SNCC’s metal prod- 
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FIG. 8 : Plot of gas mass-weighted metallicity Zi n , vs. the scaled entropy, S/Sv- Panels from left to right represent 
the cases where entropy is measured at r = 0.27?soo, r = 0.3i?soo, r = 0.5Rsoo, r = O. 8 -R 500 and r = Rsoo- Only 
clusters with Chandra coverage at each of these radii are represented in the corresponding panel. Turquoise data 
points are for CC cluster, while red points represent NCC clusters. 



FIG. 9: Plot of gas mass-weighted metallicity Zi n , vs. the scaled entropy, S/Sv- Panels from left to right represent 
the cases where entropy is measured at F g = 0.050, F g = 0.10, F g = 0.24, F g = 0.47 and F g = 0.63. Only clusters 
with Chandra coverage at each of these radii are represented in the corresponding panel. Turquoise data points are 
for CC cluster, while red points represent NCC clusters. 



FIG. 10: Distribution of the values of f ou t from 
the ensemble of Zi„ measurements. See the text in 
Section IXT1 for a description of font- 

ucts, also observ e d in e .g. iSanders fc Fabiaii (2006), 
ISimionescu et al.1 (|2009aH and iSimionescu et al ( 2010 1. 
Metal-rich gas from a small cluster, or from a galaxy, 
merging into a larger one can avoid mixing with the bulk 
ICM at large radii if its entropy is low enough to allow it 
to pierce through the outer ICM and reach the center of 
the cluster. In addition, the absense of this discrepancy 


in the outer regions of clusters could imply that met¬ 
als in the outer regions come from a source with a more 
uniform metallicity level, i.e. a source whose metallic¬ 
ity varies less from cluster to cluster. Such a homoge¬ 
neous source could be the gas that is accreted very early 
in the formation of clusters, and whose metal contribu¬ 
tion is generally referred to as pre-enrichment, which we 
further discuss in Section l4~3l Along the same lines, a 
slightly different interpretation of this finding is that, as 
time goes by, the haloes that merge later tend to have a 
wider distribution of metallicities, than those that merge 
earlier. 

The hypothesis that merger statistics is behind the 
observed metallicity distribution across clusters can be 
tested in models that combine the statistics of struc¬ 
ture formation with metal production in merging haloes, 
using a semi-an alytical approach. Such a model was 
built in lElkholvl (1201211 to test other aspects of cluster 
chemical and dynamical histories, and can be adapted to 
test whether such observations can be reproduced semi- 
analytically. 

We provide here, however, a very crude test of the 
above hypothesis. We would like to test whether the 
metallicities of clusters formed through the mergers of 
low-mass clusters from our sample have a similar distri- 
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4 DISCUSSION 


bution to that of metallicities of our high-mass clusters. 
To this end, we start by computing the metallicity Zi n re¬ 
sulting from the merger of two clusters from the sample of 
clusters in Figure[3]with M 500 < 3.5 x 10 14 M 0 . We draw 
any two clusters from the low-mass subset, and define the 
metallicity resulting from_thcir merger as the gas-mass- 
weighted average of the Zi n values of the two merging 
clusters. The gas mass used in weighting the average is 
that measured within 0.151?5oo, from the data. The re¬ 
sulting metallicity is assigned to a mass, M 500 , which is 
the sum of the masses of the merging clusters. This is 
repeated with three-four-, five- and six-cluster mergers, 
again based on our Zi n measurements in low-mass clus¬ 
ters. Our generated metallicity distribution thus comes 
from all possible mergers between 6 or less clusters from 
our low-mass clusters sample. 



m 50 „(M s ) 


FIG. 11: The width of the metallicity distribution, 
resulting from the hypothetical merger of low-mass 
clusters from our sample. The light (dark) blue re¬ 
gion represents the 95% ( 68 %) spread of generated 
metallicities. The green and orange data points 
represent our measurements in low- and high-mass 
clusters, respectively. The vertical line is plotted at 
A /500 = 3.5 x 1 0 14 Mq, and separates the above 2 
samples. 

Figure [Tl] shows the result of the above simulation, 
superimposed on our data points. We recover qualita¬ 
tively the trend in our data, whereby the width of the 
distribution of clusters’ metallicities decreases with mass. 
Beyond this toy model, a model incorporating more de¬ 
tailed statistics of structure formation, as well as a more 
sophisticated model for metal content in merging clus¬ 
ters, is required to lend support to the hypothesis linking 
our metallicity measurements to structure formation. 

4.3. Support for Pre-Enrichment? 

We argue above that the fact that metallicity disper¬ 
sion as a function of cluster mass is not observed to 
change in the outer regions, while it does vary in in¬ 
ner regions can be explained if, outside the cluster core, 
most of the metals are the result of pre-enrichment, which 
is the metallicity l evel set before cluster formation (e.g. 
Fuiita et~ahll2008t iMatsushita et al.l 120131 iWerner et al.1 
20131) . The pre-enriched gas would have an approxi¬ 
mately universal metallicity level, compared to the cold 


infalling haloes contributing metals at later times, whose 
gas metallicity values are more diverse. 

We find another clue pointing to this initial pre¬ 
enrichment metallicity level — presumed to be the same 
for all clusters — when we compare Zj n values to the 
stellar-to-gas mass fractions in clusters. iDai et all (120101 ) 
measure the stellar and baryon mass fractions of clusters 
with temperatures 1 < kTx 10 keV, using 2MASS 
data for optical measurements, and ROSAT data for X- 
ray measurements. They measure a decreasing stellar-to¬ 
gas mass ratio, r sg , as a f unction of c luster temperature, 
kTx■ Using the results of IDai et al.l (jMcj), we can thus 
estimate the stellar mass from the X-ray temperature of a 
cluster and test whether it correlates with the metallicity 
of the ICM. A correlation is expected if the population of 
stars producing the metals detected in X-ray is the same 
a s the one p r oducin g the optical luminosity of galaxies. 

IDai et akl (120101) fit the temperature-stellar-mass 
trend to a power-law and obtain the following best-fit 
relation: 


log(r Sff ) = -0.65 - 1.03 log ■ ( 14 ) 

Now, we recall that Z in = M'pJ (A Fe M™ s ), where 
M l p e and M“ s are the iron and gas masses interior to 
r = 0.15i?5oo, respectively, and Ap e is the solar iron 
abundance. We can then write 


™ M™ as A Fe 


(15) 


where M* is the stellar mass of the cluster, and 
(Mp e /M*), the average iron mass in the ICM per stellar 
mass. Or, defining Z e ff = (Mp e /M*)/A Fe , this becomes 


Z in = Z e f f r sg (kT x ) , (16) 

where Z e f f is the average iron mass in the ICM per stel¬ 
lar mass, scaled by the solar iron abundance, A Fe - The 
assumption here is that Z e ff will be the same for all 
clusters, and will not depend on kTx- 

If we assume that all the iron inside O.I 5 .R 500 has been 
produced by the stars i n the g alaxies whose luminosities 
were used by IDai et al.l (120101) to measure r sg , above, or 
that at least a population of stars of mass proportional 
to M* produced all the iron observed, then Equation ITB1 
should describe the Zi n — kTx data, once scaled by a 
suitable Z e ff. 

We perform a least-squares fit of Equation [1(1] to the 
data. The purpose of this fit is not to estimate Z e ff , but 
simply to test whether the above picture is consistent 
with our data. We assume equal errors on metallicity 
measurements so that the results are not largely biased 
by the datapoints with smaller uncertainties. The best- 
fit model plotted in Figure [12] as the dashed yellow line, 
shows that Equation [16] is an inadequate f it to the ob ¬ 
servat ions. Similar to the conclusion in iBregman et al.1 
(f2OT0l) . this suggests that the stellar population produc¬ 
ing the metals within r = O.I 5 .R 500 is unlikely to be 
related to the currently observed galaxy population in 
clusters. 

At the risk of having a model too flexible for our 
dataset, we thus introduce another parameter into Equa- 
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M 500 (M s ) kT x (keV) 

FIG. 12: Left: Zi n vs. Msoo- Right: Zi n vs. kTx- 
Red points are for NCC clusters, while turquoise 
points are for CC clusters. The yellow dashed line 
represent the Zi n — kTx best-fit based on the stellar- 
to-gas mass ratio. The purple line is a similar fit 
based on stellar-to-gas mass ratio and assuming in 
addition a baseline pre-enrichment level. 

tion Mi which is the mean metallicity in the ICM be¬ 
lieved to have been set before the current galaxy popula¬ 
tion started adding more metals. We denote this initial 
metallicity by Z pre . Our model then becomes 

Zin = Z eff r sg (kT x ) + Zpre • (17) 

Unlike Equation [TO] Equation [l7] fits the high-mass 
data better and passes near the middle of the wide distri¬ 
bution of metallicities of low-mass clusters. The best-fit 
is shown in Figure [12] as the purple line. This gives sup¬ 
port to a mass-independent initial metallicity level in the 
cores of clusters. In this picture, larger clusters would 
have most of their core metals set by pre-enrichment, 
as Z pre 3> Z e ffr sg for large kTx ■ Conversely, smaller 
clusters would have a larger contribution from metals 
associated with stars, as Z e ffr sg (kTx) increases at low 
kTx ■ The best-fit model has a pre-enrichment value of 
Z pr e ~ O.3Z 0 , and a Z e ff 2 .Zq. While the value 
of Z pre agrees with expectations, the large Z e ff value 
again points to the lack of observed galaxies, compared 
to the observed metals in clusters. The two fits above 
give similar results when repeated after excluding the 
bias-suspected clusters, i.e. clusters with f out < 0.25. 

4.4. Inner Metallicity Boost During Mergers 

Finally, we test whether the dispersion in Zi n is re¬ 
lated to dynamic activity, as measured by the centroid 
shift, (w), which is the size of the scatter of the X- 
ray centroid measured within various ap e rtures around 
the X-ray peak (e.g. iMohr et al.l 119931 : lO’Hara et al.l 
120061 IPoole et~ahl 120061) . and which is ideal for captur¬ 
ing merger activity (IPoole et al.l I2006f ) . We find weak 
evidence for a Zi n — ( w ) correlation, as we describe be¬ 
low. 

Following the prescription in IPoole et~ahl (120061b we 
calculate the position of the centroid of the X-ray emis¬ 
sion within a radius of 0.3-Rsoo, and excluding the cen¬ 
tral 30kpc. We then calculate the centroid for apertures 
that are successively smaller by 5% of 0.3i?5oo- For each 
aperture, i , we record the distance between calculated 
centroid and the X-ray peak, di. The centroid shift, ( w ), 
is then simply the standard deviation of the distances, 


di, scaled by i?soo- We compute (w) for clusters with a 
FOV that covers 0.3i?soo entirely. We find this radius to 
be a good compromise to include a large area and a good 
number of clusters. We show the computed ( w) values in 
Table lEP Figure fl3l shows a plot of Z in vs. ( w ), in the 
right panel. We also include the same plot for Z m id , in 
the left panel of the figure. Datapoints corresponding to 
low-mass clusters are in green color to distinguish them 
from the high-mass cluster datapoints in orange. 



FIG. 13: Left: Z m a vs. ( w ). Right: Zi n vs. (w). 
Green data points in either panel are for low-mass 
clusters, where M 500 < 3.5 x 10 14 Mq. Orange data 
points are for M 500 > 3.5 x 10 14 Mq clusters. 

Figure [13] suggests that there might be a correlation 
between Zi n and (w), for low-mass clusters. Pearson’s 
product-moment correlation coefficient is found to be 
rp = 0.66, with a p-value of 1.0% for the Z ln — (w) 
correlation in clusters with M 500 < 3.5 x 10 1 a Mq. 

Tins result suggests an alternative explanation to 
the Zi n dispersion in low-mass clusters, whereby merg¬ 
ers boost measured metallicity. This can be achieved 
in one of two ways. First, if mergers induce cen¬ 
tral AGN activity, then we might be observing met¬ 
als distributed into the central cluster region by the 
central active engine, as was ob serv ed and described 
in e.g . iSimiones cu et all (l2009alfQ) : iKirknatrick et al.l 
(|2009L I2011D : lO’Sullivan et al.l ( 201 11 1. This can also be 
the reason that we observe correlations between Z; ln and 
the scaled entropy, x s , as was shown in Section 13731 Sec¬ 
ond, this enhanced central metallicity might simply be 
the measurement of the central metallicity peak of an in- 
falling sub-cluster that is not completely merged with the 
main cluster, and whose emission is superimposed on the 
r < 0.151?5oo region of the main cluster. In both cases, 
these effects would need to have a stronger influence on 
metallicity in low-mass clusters compared to larger clus¬ 
ters. 

4.5. Searching for a Pre-Enrichment-Pre-Heating Link 

The motivation behind undertaking the study of metal- 
licity and entropy as far as possible from the central re¬ 
gion where the influence of the central engine and 
the effects of cooling increase — was to search for the 
signature of supernovae (SNe), which have been heating 
and chemically enriching the gas surrounding them even 
before the formation of the galaxy clusters. 

The plots of Z m id vs. x s in Section l3~3l show no hint 
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of a relation between the two quantities. However, we 
note that for a given radius, the range of values for 
x s = S/S grav can span a range as large as a factor of 
3, at large radii. One can thus envision further study 
at such large radii to be applied to a larger sample of 
clusters, which is available in current Chandra archival 
data. A stacking technique can be used to look for trends 
between metallicity and excess entropy, and to lower any 
systematics due to the high X-ray background count frac¬ 
tion at such high radii. For example, we could group 
clusters in bins of S/S grav , where S is measured at a 
large scaled radius, then extract spectra from each group 
of clusters from uniform i? 5 oo-scaled radial bins. These 
spectra could then be simultaneously fit, assuming they 
all have the same metallicity at a given scaled radial bin, 
and allowing for the temperatures to vary to match each 
cluster’s temperature. Such study could be more sensi¬ 
tive to a potential weak trend between metallicity and 
excess entropy, pointing to the effects of early supernova 
enrichment and heat injection. 

5. CONCLUSION 

In this work, we analyze a sample of 46 galaxy clusters, 
extracting chemical and dynamical measurements, in the 
hopes of obtaining clues about the history of clusters. We 
measure entropy profiles out to the largest radii where 
temperature can be measured, and provide the best-fit 
temperature and density profiles for the community to 
use. We also measure metallicity profiles, for our clus¬ 
ter sample, and present them below. The data is made 
available on an FTP sitJ^I. 

We observe a difference in the scaled iron mass between 
the centers of low-mass clusters, and the centers of high- 
mass clusters: the values of the iron content in small 
clusters are more dispersed than those in large clusters. 
We suggest two possible interpretations of this observa¬ 
tion: 

1. The lower dispersion in the larger clusters may be 
a result of the averaging of metallicities from the 
larger number of haloes that have merged to form 
them. The fact that this effect can be seen even 
in the core of clusters lends support to the idea 
that the centers of clusters continue being enriched 
by cold and metal-rich gas, originating from the 
cluster outskirts, even at low redshift. 

2 . Alternatively, there are hints that clusters can un¬ 
dergo a boost of metallicity during a merger event, 
which can contribute to the enhancement of metal¬ 
licity measured in low-mass clusters. 

We also look for a connection between the bulk metal 
content of clusters and their dynamical state, as mea¬ 
sured by the deviation of their entropy profiles from 
a self-similar profile, expected from gravitational shock 
heating, during cluster formation. We find no evidence 
of such relation in our data. More sophisticated studies 
using a larger sample would be required for such mea¬ 
surement to obtain more conclusive results. 

TYE would like to thank Michael McDonald for very 
fruitful conversations. 

10 ftp://space.mit.edu/pub/tamer/ebc2015/ 
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A ENTROPY PROFILES 


APPENDIX 

A. ENTROPY PROFILES 



r(Mpc) 


r(Mpc) 


FIG. Al: Entropy profiles of the individual clusters in the sample. The turquoise area represents the 68% confidence 
region computed at various radii, from the temperature and density profiles produced by the Monte Carlo simula- 
tions described in Sect ion 12.3.31 The gray error bars are the entropy profile measurements of the ACCEPT study, 
(ICavagnolo et al]|2009l l. See text for details. The vertical red line is drawn at r — Rsoo- 
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FIG. Al: Continued. 
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FIG. Al: Continued. 
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FIG. Al: Continued. 
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FIG. Al: Continued. 
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FIG. Bl: Metallicity profiles of observed Chandra clusters. The two vertical red lines are drawn at r 
and r = R$oo. See Section [2.3.41 for the details of the metallicity measurement method. 
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FIG. Bl: Continued. 
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Z(Solar) Z(Solar) Z(Solar) Z(Solar) 


iiizw54 


mkw3s 



0.001 0.010 0.100 0.01 0.10 1.00 
R (Mpc) R (Mpc) 


mkw8 pks0745-191 



o.ooi o.oio o.ioo 

R (Mpc) 


0.01 0.10 1.00 
R (Mpc) 


ugc3957 



zl 215 

0.5 b . . . 


ns 

o 


« 

N 


0.4 

0.3 


0.2 


0.1 


0.0 .. .. .. 

0.01 0.10 1.00 
R (Mpc) 


zl 742 



FIG. Bl: Continued. 








































































30 


C CLUSTER MASSES AND SCALES 


C. CLUSTER MASSES AND SCALES 


TABLE Cl: Cluster Masses and Scales. 


Cluster 

kTx (keV) 

7?50o(Mpc) 

M 5 oo(10 14 M 0 ) 

Mgas(W 13 M Q ) 

Abell 119 

5-79+1 

1 1 q+ 0.02 
± ‘ xo — 0.02 

4 46” 1 ” 0 ' 2 
^•^ u - 0.2 

4-50+1 

Abell 1413 

8-42+1 

i q-, +0.007 
1,0 -0.03 

7.641+ 

8 - 111 + 

Abell 1644 

c qi +0.05 
0,0 -0.09 

i i c+ 0.02 

1,1 - 0.02 

4-781+ 

5-50+1 

Abell 1651 

7-52+1 

1 OQ + 0.01 

l-*o_0.03 

6 . 8 II+ 

7-30+1 

Abell 1689 

10-4^02 

1 A5+0-005 

10.9 + 1 

12 4+ 0 - 1 
1 Z -^-0.4 

Abell 1736 

3 1Q + 0 - 05 
-0.09 

0 920 -1 ” 003 
u.yzu_o oog 

2 43 + 0 - 3 
z -^°-0.07 

2.80+1 

Abell 1795 

6 75+0-05 
°- ' D -0.05 

i oq+ 0.01 
1 - zo _ 0.02 

5-861+ 

6-19+1 

Abell 1914 

10 i +0 - 5 

1U -0.5 

1 45 + 0.004 
1 -^°-0.05 

10.7 + 1 

12.31+ 

Abell 2029 

o on+0.09 
o-^U_ 0>09 

1 44 + 0 - 06 

9.6011 

12.5+; 

Abell 2063 

3-451 + 8 

0 906 +0 03 
u.yuo _ 0 009 

9 tl +0-07 
z,ox - 0.2 

2-361+0 

Abell 2065 

6-is±g;l 

1 19+ 002 
x-xy — 0.02 

5.36+1 

5-80+1 

Abell 2142 

n-2±°;| 

1 C 7 + 0.01 

1,0 -0.06 

12 . 6+ 2 

14 5+ 0 ' 1 

Abell 2147 

4-56+1 

1 03”*” 003 
x,uo —0.009 

3.37+1 7 

3 44+0 - 1 
-0.08 

Abell 2163 

23.4:+ 

2 14+0- 07 
z - x ^-o.i 

35.1+; 

43.2+’ 

Abell 2204 

9.47+J 

1 47+0.005 
x '^ ‘ —0.06 

11 . 0+ 1 

13-6+1 

Abell 2244 

5-881+ 

1 1 O+ 0.02 

J --to_ 0.02 

5-31+1 

6-041+ 

Abell 2256 

7 7 q + 0.3 
' * '°-0.3 

1 q 7 +0.03 
xo -0.03 

8 - 11+1 

9.541+ 

Abell 2319 

10 I "*" 0,3 
lu ‘- 0.2 

4-581+0 

12 .5 + 3 

15.61+ 

Abell 2657 

3.821+ 

0.902l+° 3 4 

2.28+1 4 

2 - 081+7 

Abell 2734 

4-28+2 

0-982++7 

3.00+1 5 

3 02”*” 0 ' 1 
°- uz -0.08 

Abell 3112 

5.01+- 

1 09”*” 002 
1 ,uy _ 0.02 

4 i9 +°- 2 
^• iy - 0.2 

4.81+1 

Abell 3158 

5-2o±8:8S 

1 in - * -002 
x - xu -o.oi 

4.23+1 

4-531+g 

Abell 3376 

4-581+ 

0 964”*" 0,03 

U.»U^_Q QQ 2 

2.801+ 3 

9 49+ 0 08 

Abell 3391 

5 41+ 0 - 4 
O.^l_o 2 

1 in - * -002 
xxu - 0.02 

4.19+; 3 

4-28+1 

Abell 3571 

7 iyq + 0.3 
' • '°-0.3 

1 qc+ 0.01 

T.OO-g 03 

7.60+1 

8-431+ 

Abell 3667 

6-60++ 

1 qo+ 0-01 

± -°^_0.03 

7 qq+ 0.2 

< ••J'J-0.5 

9.381+ 

Abell 3822 

5.361+ 

1 1 c+ 0.02 
1 - lo - 0.02 

4.82+1 

5-571+ 

Abell 3827 

7-891+ 

1 qq+0-01 

t.oo_o 04 

7.60+; 3 

8.461+ 

Abell 3921 

5-931+ 

1 1 q+ 0.02 

1 - i 8 - 0.02 

5.39+1 

6-15+1 

Abell 399 

6-471+ 

1 - 261+1 

6.33+1 

7 43 + 0- 1 

Abell 400 

9 1 5+0-05 
z,lo -0 04 

0.6781+1 

0.955++ 

Q 799 +O O 3 
u. 1 ^»_o.07 

Abell 4038 

o -i q+0.05 
0.03 

0.839l+°3 8 

1.82+1 5 

1 7-1 +0.04 
x ‘' x —0.09 

Abell 4059 

4-341+2 

0 948” 1 " 0,03 
u.» i ±o_ 0 004 

2 66”*” 0 ' 2 
zoo -0 03 

2 4i+ 0 05 
z '^ x — 0.07 

Abell 478 

7-651+ 

1 qn+0.006 

L.oy_ 0 04 

8 . 801 + 

11 3 +04 
xx -°- 0.2 

Abell 539 

2 59+ 0 - 04 
^•oy_ 0 .04 

0-8031+1 

1-59+1 5 

1 63”*” 0 ' 1 

T.UO-0 03 

Abell 644 

8.491+ 

1-301+1 

7.06l+ 5 

6 - 861 + 

Abell 754 

11 . 81 + 

1 qq+0.05 
T.OO-0 04 

7-441+ 

5-381+ 

Abell S 405 

4-621+ 

0.985++7 

3-031 + 7 

2-85+1 9 

Hydra A 

q 7 C+0.04 
' <J -0.03 

1.021+08 

3.36++ 8 

4.201 + 7 

Zw III 54 

9 95+0-06 
^•^-0.06 

0.7181+1 

1 14+0-07 
1 - 1 ^- 0.2 

1 05+0-03 

i.UO_o 09 

MKW 3s 

3 44 + 0 09 
°*^-0.03 

0.909l+° 3 4 

9 qc+0.03 

z.oo_ 0 2 

2.43+1 3 

MKW 8 

2.50+1 

0 715+ 0 02 
u. t J-o_ 0 04 

1 1 q+0-07 
x - xo -0.2 

0 924 -1 "® 05 
u.yz ^_ 0 07 

PKS 0745-191 

6-761+ 

1 qc+0.01 

l.OO_ 0 04 

8.09+1 

ii-il+ 

UGC 3957 

2.34+1 

0 697~*~ 0 ' 02 
uoy * -0.04 

1-051 + 7 

0-864+;“ 

ZwCl 1215+0400 

7 ^7+0-3 
‘ - 0 ' -0.3 

1 97+0.01 

± - z ‘ -0.03 

6-57++ 

6-78+1 

ZwCl 1742+3306 

4-461+ 

0 957 ~*~ 003 

u.yo / —0.006 

2.811+ 6 

2.611+s 
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D. DENSITY AND TEMPERATURE PROFILES BEST-FIT PARAMETERS 


TABLE Dl: Best-fit parameters of the electron density radial profiles. 


Cluster 

n 0 

a 

P 

r c 

r s 

e 

Abell 119 

0.0005406 

0.627 

5.0 

2.959 

0.289 

1.227 

Abell 1413 

0.04096 

0.0 

0.3753 

0.0217 

0.3968 

2.13 

Abell 1644 

0.04399 

0.9245 

0.3237 

0.004989 

2.17 

5.0 

Abell 1651 

0.009126 

0.6828 

0.3982 

0.08472 

0.2506 

1.249 

Abell 1689 

0.04991 

0.0 

0.399 

0.0306 

0.3455 

1.99 


Note. — Table [Dll is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. The n e (r) model is shown in Equation [l] 


TABLE D2: Best-fit parameters of the temperature radial profiles. 


Cluster 

n 

a 

b 

c 

a cool 

r cool 

tmin 

to 

Abell 119 

0.2912 

-0.3165 

4.758 

1.011 

12.51 

0.08264 

19.22 

8.419 

Abell 1413 

0.01754 

0.4569 

6.889 

-0.7317 

-0.7636 

0.1653 

1.329 

7.881 

Abell 1644 

0.4703 

-0.2443 

7.74 

0.5485 

3.737 

0.02752 

4.614 

6.027 

Abell 1651 

0.04094 

0.04568 

4.947 

-0.0793 

1.9 

0.05 

6.693 

6.858 

Abell 1689 

1.046 

-0.002614 

7.542 

3.02 

8.231 

0.07014 

10.04 

11.45 


Note. — Table [D2l is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. The kT(r) model is shown in Equation [5] 


E. MORPHOLOGICAL PARAMETERS 


TABLE El: Morphological Parameters and Entropy near the Center. 


Cluster 

css a 

•S , 40 b (keV cm 2 ) 

(w) c /10~ 3 

Abell 119 

Abell 1413 

Abell 1644 

Abell 1651 

Abell 1689 

2.17± 0.037 

9.91± 0.079 

6.13± 0.057 

7.71± 0.16 

12.3± 0.064 

550+®l; 

95.4+l'f 

140+S 

108.11? 

1.12 

0.607 

(4.60) 

3.15 

1.16 


Note. — Table [Ell is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. Values of ( w ) in parentheses correspond to measurements in observations where the FOV does not fully 
cover the region r < O. 3 R 50 O 5 an d are excluded from analysis involving ( w ), but shown here. 


a Surface brightness concentration. 
b Entropy at r =40 kpc. 
c Centroid shift. 
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F GLOBAL METALLICITY MEASURES 


F. GLOBAL METALLICITY MEASURES 


TABLE FI: Global metallicity measures. 


Cluster 

Zmid(Z q) 

Zin(Z 0 ) 

Abell 

119 

0 260 +0 ' 036 
u.zou _ 0 035 

fl QQ7+0-047 
0.001 _ 0.045 

Abell 

1413 

fl oc:o+0.032 

u.zuo_0 o 3 2 

O Q79+0-027 

u.o/z_o 026 

Abell 

1644 

0 226~*~ 0 ' 029 
u.zzo_ 0 02g 

0 30Q+O O 28 

U.OOU _ 0 026 

Abell 

1651 

n Q4S+0-081 

U.o4o_ 0 0 y 6 

0 3«3+0-069 

u.ooo _ 0 064 

Abell 

1689 

0 304+ 0 - 032 

U.OLE±_o 031 

O QK9+0-025 

u.ooz _ 0 024 

Abell 

1736 

n 9 qq+0-036 

u.zoo_ 0 034 

0 306 -1-0 ' 046 

U.OUO_o 044 

Abell 

1795 

0 240~*~ 0 ' 023 

U.z^tu _ 0 023 

0 324~*~°' 013 
u.oz "±_ 0 013 

Abell 

1914 

0 314~*~ 0 ' 095 
u - OJ - -0.091 

0 404~*~°' 064 

u.^±u^_o 061 

Abell 

2029 

n oqo+0.024 

u.zoo _ 0 023 

0 421~*~°' 015 

U.4:Z±_o 015 

Abell 

2063 

f) 1 Q7+0-024 

U.18/_0 022 

n Q 97 + 0-022 

u.oz / _ 0 021 

Abell 

2065 

H 991 +0.043 
u.zzL _ 0 Q40 

0 316+0- 041 

U.OLU _0 039 

Abell 

2142 

0 366~*~ 0 ' 037 

U.oOO—o_035 

n Q 90 + 0.023 
u.ozo _ 0 022 

Abell 

2147 

0-290t“;°“ 

0 269~*~°' 039 

u-zu»_o .037 

Abell 

2163 

0.476l°;i“ 

0 345 ~*~°' 077 
u.o^o _ 0 075 

Abell 

2204 

0 3 Q 7 +O O 43 

u.oy / —0,042 

0.39lt“;°2l 

Abell 

2244 

0.2011^25 

0.298±g;“S 

Abell 

2256 

0.289±H?8 

0.389t“°®® 

Abell 

2319 

0.286±°;°“ 

O Q77+0-049 

1 -0.047 

Abell 

2657 

fl 981 +0.043 
U.zol_o 040 

0 460+ 0 055 

U. z ±DU_o 051 

Abell 

2734 

n 997 + 0.056 
u.z^ / —0.052 

O 997+0.056 
u.zzi / —0.052 

Abell 

3112 

fl 1 OO+0.034 

u.ioz _ 0 032 

0 3Q4 + 0,029 
u.oy^:_o 027 

Abell 

3158 

0.389tg;gS 

0 42f)~*~ 0-029 
u -^ZU_o 027 

Abell 

3376 

H 91 O+0.042 
u.oiz_o 040 

0 22‘5 "*~ 0-044 

u.zzo _ 0 042 

Abell 

3391 

0 99Q+0-069 

U.zzo _ 0 064 

n 990 + 0 .O 69 
u.zzo _ 0 064 

Abell 

3571 

0 1 Q 1 +0-034 
U.iyi_o 032 

0 3Q2 + 0 - 024 

u.oyz _ 0 023 

Abell 

3667 

0 309+0- 010 
u.9uy _ 0 010 

0 343 + 0-012 
U.040— 0.011 

Abell 

3822 

0 258"*”°‘ 067 

u.zoo _ 0 063 

0 9^Q+0-067 

u.zoo_o 063 

Abell 

3827 

0 9R3+0-052 

U.ZOO _ 0 049 

0-368±S 

Abell 

3921 

n 097 + 0.051 

u.oz / —0,049 

0.308t“;°« 

Abell 

399 

n 9QO+0.047 
U.zoo_o 045 

0.307±g;°S 

Abell 

400 

0 394+0 040 

u.o »^_0 036 

O c:/ 1Q+0-059 
U.548_o 053 

Abell 

4038 

0 q^^+O-025 
u.ooo _ 0 023 

0 440+ 0,020 
U.^U_o 019 

Abell 

4059 

0 250" 1 ” 0 ' 017 
U.ZOU_Q 016 

0-4661°;°“ 

Abell 

478 

fl 9^7+0.025 

U.zo /_ 0 024 

0.328±g;°i« 

Abell 

539 

0 91 9+0.023 

U.Z1Z _ 0 021 

0 993+0-024 

u.zzo _ 0 022 

Abell 

644 

n 090 + 0.044 

u.ozo _ 0 042 

O 4QO+0.038 

U.4oo_o 036 

Abell 

754 

0 3 Q 2+0 0 58 

u.ouz —0057 

0 369~*~°' 036 
u.ou ^_ 0 035 

Abell 

S 405 

0 264 +0 ‘ 083 
u.zo4t _ 0 077 

0 264~*~°' 083 
u.ZD^_o 077 

Hydra A 

fl 1 4 7+0.0095 
-0.0094 

n 90 C+ 0.0055 
u.zou _ 0 .0054 

Zw III 54 

0 258"*” 0 ' 029 

U. zoo—o_o26 

0 370 -1-0 ' 037 

U.o / U—0.034 

MKW 3s 

n 090 + 0.025 
u.ozO_o 024 

0 395 ~*~°' 017 

u.ot 7 u_ 0 017 

MKW 8 

0 241 +0 ' 045 
u -^ x -o.039 

0 466" 1 " 0 ' 082 

U.-4:00 _o 075 

PKS 0745-191 

0 206 + °' 045 
u.zuo _ 0 043 

0 369~*~°' 028 
u.ouy_o 026 

UGC 3957 

0 213+0-051 

U.Z10_o 044 

0 426~*~°' 049 

u.^zo _ 0 043 

ZwCl 1215+0400 

pi 91 r+0.060 

U.O1O-0 056 

fl QK7+0.069 
U.oO 1 —0.065 

ZwCl 1742+3306 

0 314+0- 045 
u.O-L^-o 042 

0 447+0 - 026 
_o.o25 





